\relax 
\providecommand\HyperFirstAtBeginDocument{\AtBeginDocument}
\HyperFirstAtBeginDocument{\ifx\hyper@anchor\@undefined
\global\let\oldcontentsline\contentsline
\gdef\contentsline#1#2#3#4{\oldcontentsline{#1}{#2}{#3}}
\global\let\oldnewlabel\newlabel
\gdef\newlabel#1#2{\newlabelxx{#1}#2}
\gdef\newlabelxx#1#2#3#4#5#6{\oldnewlabel{#1}{{#2}{#3}}}
\AtEndDocument{\ifx\hyper@anchor\@undefined
\let\contentsline\oldcontentsline
\let\newlabel\oldnewlabel
\fi}
\fi}
\global\let\hyper@last\relax 
\gdef\HyperFirstAtBeginDocument#1{#1}
\providecommand\HyField@AuxAddToFields[1]{}
\citation{bednorz1986possible,anderson1987resonating}
\citation{halperin1993theory}
\citation{anderson1973resonating,elser1989nuclear,zhou2017quantum,jia2020mott,ji2021rare}
\citation{kim1996observation,giamarchi2003quantum,auslaender2005spin,jompol2009probing}
\citation{imada1998metal,lee2006doping}
\citation{tomonaga1950remarks,luttinger1963exactly,gogolin2004bosonization}
\citation{bethe1931theorie,wang2015off}
\citation{anderson1987resonating,maier2000d,maier2005systematic,gull2013superconductivity,su1998nonexistence,qin2020absence}
\citation{baxter2016exactly}
\citation{Kitaev2006}
\citation{Chen2018}
\citation{ezawa2018exact}
\citation{Miao2019}
\newlabel{FirstPage}{{}{1}{\relax }{section*.1}{}}
\@writefile{toc}{\contentsline {title}{Phase transition in BCS-Hubbard model}{1}{section*.2}}
\@writefile{toc}{\contentsline {abstract}{Abstract}{1}{section*.1}}
\@writefile{toc}{\contentsline {section}{\numberline {I}Introduction}{1}{section*.3}}
\@writefile{toc}{\contentsline {section}{\numberline {II}Model and Method}{1}{section*.4}}
\newlabel{sec:model_method}{{II}{1}{\relax }{section*.4}{}}
\@writefile{toc}{\contentsline {subsection}{\numberline {A}Field BCS-Hubbard model}{1}{section*.5}}
\newlabel{sec:model}{{II\tmspace  +\thinmuskip {.1667em}A}{1}{\relax }{section*.5}{}}
\newlabel{eqHamiltonian}{{1}{1}{\relax }{equation.2.1}{}}
\newlabel{eqHoppingHamiltonian2}{{2a}{1}{\relax }{equation.2.2alph1}{}}
\newlabel{eqPairingHamiltonian2}{{2b}{1}{\relax }{equation.2.2alph2}{}}
\newlabel{eqUHamiltonian2}{{2c}{1}{\relax }{equation.2.2alph3}{}}
\newlabel{eqMagHamiltonian2}{{2d}{1}{\relax }{equation.2.2alph4}{}}
\citation{Miao2019}
\citation{Xiao2009}
\newlabel{eqHoppingPairingHamiltonian}{{4}{2}{\relax }{equation.2.4}{}}
\newlabel{eqHamiltonianU}{{5}{2}{\relax }{equation.2.5}{}}
\newlabel{eqMagHamiltonian}{{6}{2}{\relax }{equation.2.6}{}}
\newlabel{eqCommuativeHamiltonian}{{7}{2}{\relax }{equation.2.7}{}}
\newlabel{eqConservedQuantities}{{8}{2}{\relax }{equation.2.8}{}}
\newlabel{eqDefSubspace}{{9}{2}{\relax }{equation.2.9}{}}
\newlabel{eqRefermionization}{{10}{2}{\relax }{equation.2.10}{}}
\newlabel{eqMain}{{11}{2}{\relax }{equation.2.11}{}}
\@writefile{toc}{\contentsline {section}{\numberline {III}Results}{2}{section*.6}}
\newlabel{exact_solution}{{III}{2}{\relax }{section*.6}{}}
\@writefile{toc}{\contentsline {subsection}{\numberline {A}Analytical method}{2}{section*.7}}
\newlabel{eqGeneralH}{{12}{2}{\relax }{equation.3.12}{}}
\newlabel{eqGeneralMatrix}{{13}{2}{\relax }{equation.3.13}{}}
\citation{ezawa2018exact}
\citation{ezawa2018exact}
\citation{Kitaev2006}
\citation{lieb2004flux}
\newlabel{eqXN}{{14}{3}{\relax }{equation.3.14}{}}
\newlabel{eqPicdfermions}{{17}{3}{\relax }{equation.3.17}{}}
\newlabel{eqGS}{{18}{3}{\relax }{equation.3.18}{}}
\newlabel{eqExcitationSpectrum}{{19}{3}{\relax }{equation.3.19}{}}
\@writefile{toc}{\contentsline {subsection}{\numberline {B}Phase diagram}{3}{section*.8}}
\newlabel{eqDefSectors}{{21}{3}{\relax }{equation.3.21}{}}
\newlabel{eqGap}{{22}{3}{\relax }{equation.3.22}{}}
\newlabel{eqOrderParameters}{{23}{3}{\relax }{equation.3.23}{}}
\@writefile{lof}{\contentsline {figure}{\numberline {1}{\ignorespaces The phase diagram in large scale. Here we set $t = 1.0,\Delta = 0.0, N_x = N_y = 40$, $U$ is the Hubbard interaction and $h$ is the magnetic field applied along the $y$-axis. Three different phase is distinguished by the order parameter is $D^{(+)}$ defined in (\ref  {eqOrderParameters}).}}{4}{figure.1}}
\newlabel{figPhaseDiagram}{{1}{4}{The phase diagram in large scale. Here we set $t = 1.0,\Delta = 0.0, N_x = N_y = 40$, $U$ is the Hubbard interaction and $h$ is the magnetic field applied along the $y$-axis. Three different phase is distinguished by the order parameter is $D^{(+)}$ defined in (\ref {eqOrderParameters})}{figure.1}{}}
\@writefile{toc}{\contentsline {section}{\numberline {IV}Numerical analysis}{4}{section*.9}}
\newlabel{MC_scaling}{{IV}{4}{\relax }{section*.9}{}}
\@writefile{toc}{\contentsline {subsection}{\numberline {A}scaling analysis}{4}{section*.10}}
\newlabel{SA}{{IV\tmspace  +\thinmuskip {.1667em}A}{4}{\relax }{section*.10}{}}
\@writefile{lof}{\contentsline {figure}{\numberline {2}{\ignorespaces The scaling analysis of the critical points $h_c$ for $t = 1.0, \Delta = 0.0, U = 18.0$ , the $h_c$ quickly converge to the value in thermodynamic limit when the size of lattice increases. Here, $h_{c}$ stands for the strength of critical magnetic field of the three phase boundaries.}}{4}{figure.2}}
\newlabel{figphaseregion}{{2}{4}{The scaling analysis of the critical points $h_c$ for $t = 1.0, \Delta = 0.0, U = 18.0$ , the $h_c$ quickly converge to the value in thermodynamic limit when the size of lattice increases. Here, $h_{c}$ stands for the strength of critical magnetic field of the three phase boundaries}{figure.2}{}}
\@writefile{toc}{\contentsline {subsection}{\numberline {B}Thermal properties}{5}{section*.11}}
\newlabel{TP}{{IV\tmspace  +\thinmuskip {.1667em}B}{5}{\relax }{section*.11}{}}
\newlabel{eqmy}{{24}{5}{\relax }{equation.4.24}{}}
\@writefile{lof}{\contentsline {figure}{\numberline {3}{\ignorespaces $\delimiter "426830A M_{y}\delimiter "526930B _T$ as a function of $h$ at different temperature with $\beta =10,100$ and $1000$ in a $4\times 4$ square lattice, respectively. (a)Relationship between $h$ and $\delimiter "426830A M_{y}\delimiter "526930B _T$ for small $U$ with $U=2$. (b) Relationship between $h$ and $\delimiter "426830A M_{y}\delimiter "526930B _T$ for large $U$ with $U=18$.}}{5}{figure.3}}
\newlabel{figmy}{{3}{5}{$\langle M_{y}\rangle _T$ as a function of $h$ at different temperature with $\beta =10,100$ and $1000$ in a $4\times 4$ square lattice, respectively. (a)Relationship between $h$ and $\langle M_{y}\rangle _T$ for small $U$ with $U=2$. (b) Relationship between $h$ and $\langle M_{y}\rangle _T$ for large $U$ with $U=18$}{figure.3}{}}
\@writefile{toc}{\contentsline {section}{\numberline {V}Conclusion and Summary}{5}{section*.12}}
\newlabel{sec:summary}{{V}{5}{\relax }{section*.12}{}}
\bibdata{mainNotes,bcs_hubbard}
\bibcite{bednorz1986possible}{{1}{1986}{{Bednorz\ and\ M{\"u}ller}}{{}}}
\bibcite{anderson1987resonating}{{2}{1987}{{Anderson}}{{}}}
\bibcite{halperin1993theory}{{3}{1993}{{Halperin\ \emph  {et~al.}}}{{Halperin, Lee,\ and\ Read}}}
\bibcite{anderson1973resonating}{{4}{1973}{{Anderson}}{{}}}
\bibcite{elser1989nuclear}{{5}{1989}{{Elser}}{{}}}
\bibcite{zhou2017quantum}{{6}{2017}{{Zhou\ \emph  {et~al.}}}{{Zhou, Kanoda,\ and\ Ng}}}
\bibcite{jia2020mott}{{7}{2020}{{Jia\ \emph  {et~al.}}}{{Jia, Gong, Liu, Zhao, Dong, Dai, Li, Lei, Yu, Zhang \emph  {et~al.}}}}
\bibcite{ji2021rare}{{8}{2021}{{Ji\ \emph  {et~al.}}}{{Ji, Sun, Cai, Wang, Sun, Ren, Zhang, Jin,\ and\ Zhang}}}
\bibcite{kim1996observation}{{9}{1996}{{Kim\ \emph  {et~al.}}}{{Kim, Matsuura, Shen, Motoyama, Eisaki, Uchida, Tohyama,\ and\ Maekawa}}}
\bibcite{giamarchi2003quantum}{{10}{2003}{{Giamarchi}}{{}}}
\bibcite{auslaender2005spin}{{11}{2005}{{Auslaender\ \emph  {et~al.}}}{{Auslaender, Steinberg, Yacoby, Tserkovnyak, Halperin, Baldwin, Pfeiffer,\ and\ West}}}
\bibcite{jompol2009probing}{{12}{2009}{{Jompol\ \emph  {et~al.}}}{{Jompol, Ford, Griffiths, Farrer, Jones, Anderson, Ritchie, Silk,\ and\ Schofield}}}
\bibcite{imada1998metal}{{13}{1998}{{Imada\ \emph  {et~al.}}}{{Imada, Fujimori,\ and\ Tokura}}}
\bibcite{lee2006doping}{{14}{2006}{{Lee\ \emph  {et~al.}}}{{Lee, Nagaosa,\ and\ Wen}}}
\bibcite{tomonaga1950remarks}{{15}{1950}{{Tomonaga}}{{}}}
\bibcite{luttinger1963exactly}{{16}{1963}{{Luttinger}}{{}}}
\bibcite{gogolin2004bosonization}{{17}{2004}{{Gogolin\ \emph  {et~al.}}}{{Gogolin, Nersesyan,\ and\ Tsvelik}}}
\bibcite{bethe1931theorie}{{18}{1931}{{Bethe}}{{}}}
\bibcite{wang2015off}{{19}{2015}{{Wang\ \emph  {et~al.}}}{{Wang, Yang, Cao,\ and\ Shi}}}
\bibcite{maier2000d}{{20}{2000}{{Maier\ \emph  {et~al.}}}{{Maier, Jarrell, Pruschke,\ and\ Keller}}}
\bibcite{maier2005systematic}{{21}{2005}{{Maier\ \emph  {et~al.}}}{{Maier, Jarrell, Schulthess, Kent,\ and\ White}}}
\bibcite{gull2013superconductivity}{{22}{2013}{{Gull\ \emph  {et~al.}}}{{Gull, Parcollet,\ and\ Millis}}}
\bibcite{su1998nonexistence}{{23}{1998}{{Su\ and\ Suzuki}}{{}}}
\bibcite{qin2020absence}{{24}{2020}{{Qin\ \emph  {et~al.}}}{{Qin, Chung, Shi, Vitali, Hubig, Schollw{\"o}ck, White, Zhang \emph  {et~al.}}}}
\bibcite{baxter2016exactly}{{25}{2016}{{Baxter}}{{}}}
\bibcite{Kitaev2006}{{26}{2006}{{Kitaev}}{{}}}
\bibcite{Chen2018}{{27}{2018}{{Chen\ \emph  {et~al.}}}{{Chen, Li,\ and\ Ng}}}
\bibcite{ezawa2018exact}{{28}{2018}{{Ezawa}}{{}}}
\bibcite{Miao2019}{{29}{2019}{{Miao\ \emph  {et~al.}}}{{Miao, Xu, Zhang,\ and\ Zhang}}}
\bibcite{Xiao2009}{{30}{}{{wen Xiao}}{{}}}
\bibcite{lieb2004flux}{{31}{2004}{{Lieb}}{{}}}
\@writefile{toc}{\contentsline {section}{\numberline {VI}Acknowledgments}{6}{section*.13}}
\@writefile{toc}{\contentsline {section}{\numberline {}References}{6}{section*.14}}
\newlabel{LastBibItem}{{31}{6}{\relax }{section*.14}{}}
\@writefile{lof}{\contentsline {figure}{\numberline {4}{\ignorespaces The results of the numeric ergodic calculation for the $4 \times 4$ lattice with $t = 1.0, \Delta = 0$. The order parameter here is $D^{(+)} = \DOTSB \sum@ \slimits@ _l D_l$.}}{7}{figure.4}}
\newlabel{figCombine}{{4}{7}{The results of the numeric ergodic calculation for the $4 \times 4$ lattice with $t = 1.0, \Delta = 0$. The order parameter here is $D^{(+)} = \sum _l D_l$}{figure.4}{}}
\@writefile{toc}{\appendix }
\@writefile{toc}{\contentsline {section}{\numberline {A}Numerical results}{7}{section*.15}}
\newlabel{appendixc}{{A}{7}{\relax }{section*.15}{}}
\@writefile{toc}{\contentsline {subsection}{\numberline {1}Configuration}{7}{section*.16}}
\@writefile{toc}{\contentsline {subsection}{\numberline {2}Pseudo spin gap and pseudo charge gap}{7}{section*.17}}
\@writefile{lof}{\contentsline {figure}{\numberline {5}{\ignorespaces The results of the pseudo spin gap and pseudo charge gap for the $4 \times 4$ lattice with $t = 1.0, \Delta = 0$. The size effect is remarkable as the system size is small. (a), the pseudo spin gap. (b), the pseudo charge gap. As presented in (a) and (b) in some regions, the pseudo spin gap is close while the pseudo charge gap is open.}}{7}{figure.5}}
\newlabel{fig_sp_and_ch}{{5}{7}{The results of the pseudo spin gap and pseudo charge gap for the $4 \times 4$ lattice with $t = 1.0, \Delta = 0$. The size effect is remarkable as the system size is small. (a), the pseudo spin gap. (b), the pseudo charge gap. As presented in (a) and (b) in some regions, the pseudo spin gap is close while the pseudo charge gap is open}{figure.5}{}}
\bibstyle{apsrev4-1}
\citation{REVTEX41Control}
\citation{apsrev41Control}
\newlabel{LastPage}{{}{8}{}{}{}}
